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We report the observation of the flavor-changing neutral current process b — > d'y using a sample of 
386 x 10 6 B meson pairs accumulated by the Belle detector at the KEKB e + e collider. We measure 
branching fractions for the exclusive modes B~ — > p~7, B° — » p°7 and B° — > wy. Assuming that 
these three modes are related by isospin, we find B(B — > (p,o;)7) = (1.32 ^'^(stat.) ioJg(syst.)) x 
10 -6 with a significance of 5.1cr. This result is used to determine the ratio of CKM matrix elements 
\V td /V ts \ to be 0.199 ^;°25 ( ex P-)taoi5( th eo.). 

PACS numbers: 11.30.Hv, 13.40.Hq, 14.65.Fy, 14.40.Nd 



The b — > d'y process, which proceeds via a loop diagram 
(Fig. HJa)) in the Standard Model (SM), is suppressed 
with respect to b — > S7 by the Cabibbo-Kobayashi- 
Maskawa (CKM) factor 1] \V t d/V ts \ 2 ~ 0.04, with large 
uncertainty due to the lack of precise knowledge of \Vtd\- 
The exclusive modes B —* p7 and B° — * ujj are presum- 
ably the easiest modes to search for; no evidence for the 
decays has been previously reported |2|,|3|. The predicted 
branching fractions are (0.9-2.7) x 10 -6 0, Q based on 
the measured rate for the b — > process B — ► K*j and 
the |Vtd/Vts| 2 factor with corrections due to form factors, 
SU(3) breaking effects, and, for the B~ decay, inclusion 
of an annihilation diagram (Fig. ^b)). Measurement of 
these exclusive branching fractions allows one to deter- 
mine the value of |T4d/V* s | in the context of the SM and 
to search for physics beyond the SM In this Letter, 
we report the observation of the b — > d'y process using a 
sample of (386 ± 5) x 10 6 B meson pairs accumulated at 
the T(4S) resonance. With a larger data sample and an 
improved analysis procedure, the results supersede those 
of our previous publication [2| . 

The data are produced in e + e~ annihilation at the 
KEKB energy-asymmetric (3.5 on 8 GeV) collider and 
collected with the Belle detector [8j, which includes a 



(a) loop diagram (b) annihilation diagram 




FIG. 1: (a) Loop diagram for b — » d'y and (b) annihilation 
diagram, which contributes only to B~ — » p~7- 



silicon vertex detector (SVD), a central drift chamber 
(CDC), aerogel threshold Cherenkov counters (ACC), 
time-of-fiight (TOF) scintillation counters, and an elec- 
tromagnetic calorimeter (ECL) of CsI(Tl) crystals lo- 
cated inside a 1.5 T superconducting solenoid coil. 

We reconstruct three signal modes, B~ — * p~*y, B° — > 
p°7 and B° — > ujj, and two control samples, B~ — ► K*~j 
and B° — > K*°j. Charge conjugate modes are implicitly 
included throughout this Letter. The following decay 
modes are used to reconstruct the intermediate states: 

p~ — > 7T~7T , p° — > 7T + 7T~, U) — > 7T + 7r _ 7r , K*~ — » K~1T , 

K*° — > K~ir + , and 7r° — * 77. 
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Photon candidates are reconstructed from ECL en- 
ergy clusters with a photon-like shape and no asso- 
ciated charged track. A photon in the barrel ECL 
(33° < 8 1 < 128° in the laboratory frame polar an- 
gle) with a center-of-mass (cm.) energy in the range 
1.8 GeV < E 1 < 3.4 GeV is selected as the primary pho- 
ton candidate. To suppress backgrounds from ir° /r] — ► 77 
decays, we apply a veto algorithm based on likelihoods 
to be and not to be a tt° /t). The likelihoods are calcu- 
lated for every combination of the primary photon and 
another photon in the event using the energy of the other 
photon and the invariant mass of the pair. We also reject 
the primary photon candidate if the ratio of the energy 
in the central 3x3 ECL cells to that in the central 5x5 
cells is less than 0.95. 

Neutral pions are formed from photon pairs with in- 
variant masses within ±16 MeV/c 2 (~3cr) of the ir° mass. 
The photon momenta are then recalculated with a tt° 
mass constraint. We require the energy of each photon 
to be greater than 50 (100) MeV inside (outside) the bar- 
rel ECL. We also require the cosine of the angle between 
the two photons in the laboratory frame to be greater 
than 0.7; this requirement suppresses the copious combi- 
natorial background with momenta below 0.6 GeV/c. 

Charged pions and kaons are selected from tracks in the 
CDC and SVD. Each track is required to have a trans- 
verse momentum greater than 100 MeV/c and a distance 
of closest approach to the interaction point of less than 
0.5 cm in radius and ±3.0 cm along the z-axis, which is 
parallel to the positron beam. We do not use a track to 
form the signal candidate if, when it is combined with 
an oppositely charged track, the resulting pair has an 
invariant mass within ±30 MeV/c 2 of the Kg mass and 
a displaced vertex that is consistent with that of a Kg. 
We determine pion (£„) and kaon (£k) likelihoods from 
ACC, CDC and TOF information and form a likelihood 
ratio £ 7r /(£ 7r + Lk) to separate pions from kaons. The 
criteria for pions have efficiencies of 83%, 81% and 91% 
for p~ , p and u, respectively; the corresponding kaon 
misidentffication rates are 5.8%, 6.3% and 8.4%. For K* 
candidates, we select kaons with an efficiency of 90%. 

Invariant masses for the p, to and K* candidates 
are required to be within windows of ±150, ±30 and 
±75 MeV/c 2 , respectively, around their nominal values. 

Candidate B mesons are reconstructed by combining 
a p or oj candidate with the primary photon and calcu- 
lating two variables: the beam-energy constrained mass 
M bc = v^bcam/c 2 ) 2 - IPb/ c I 2 ' and the energy differ- 
ence AE = E* B - ££ eam . Here, p B and Eg arc the 
cm. momentum and energy of the B candidate, and 
-^beam ' s the c - m - beam energy. To improve resolution, 
the magnitude of the photon momentum is replaced by 
(^beam ~ ^*p/uj)l c when the momentum p B is calculated. 

To optimize the event selection, we study Monte Carlo 
(MC) events in a signal box defined as 5.273 GeV/c 2 < 
M hc < 5.285 GeV/c 2 and -0.10 GeV < AE < 0.08 GeV. 



We choose selection criteria to maximize Ns/^/Nb, 
where TVs an d Nb are the expected signal and the sum 
of the background yields. 

The dominant background arises from continuum 
events (e + e~ — > qq{^), q = u,d,s,c), where a random 
combination of a p or u> candidate with a photon forms 
a B candidate. We suppress this background using the 
following quantities: (1) T, a Fisher discriminant con- 
structed from 16 modified Fox- Wolfram moments [fll [HI] 
and the scalar sum of the transverse momenta of all 
charged tracks and photons. (2) cos9* B , where 9* B is the 
cm. polar angle of the B candidate direction: true B 
mesons follow a 1 — cos 2 9* B distribution, while candi- 
dates in the continuum background are almost uniformly 
distributed. (3) Az, the separation along the z-axis be- 
tween the decay vertex of the candidate B meson and the 
fitted vertex of the remaining tracks in the event. Dis- 
crimination is provided due to the displacement of the 
signal B decay vertex from the other B, as tracks from 
continuum events typically have a common vertex. For 
each of the quantities J 7 , cos9 B and Az, we construct 
likelihood distributions for signal and continuum events. 
The cos 9* B and signal Az distributions are determined 
from MC samples; the continuum Az distribution is de- 
termined from the data sideband 5.20 GeV/c 2 < Mb c < 
5.24 GeV/c 2 , -0.1 GeV < AE < 0.5 GeV. 

We form product likelihoods C s and C c for signal and 
continuum background, respectively, from the likelihood 
distributions for T, cos9* B and (where available) Az. In 
addition, we use a tagging quality variable r that indi- 
cates the level of confidence in the B-flavor determina- 
tion as described in Ref. I n the (r, TZ) plane de- 
fined by the tagging quality r and the likelihood ratio 
TZ = C s / (C s + C c ), signal tends to populate the edges at 
r = 1 and TZ = 1, while continuum preferentially popu- 
lates the edges at r — and TZ — 0. We divide the events 
into six bins of r (two bins between and 0.5, and four 
between 0.5 and 1) and determine the minimum TZ re- 
quirement for each bin. In the p~"f mode, we also assign 
events to the bin < r < 0.25 if the tagging-side flavor 
is the same as the signal-side. The signal efficiency is 
^40%, and ~95% of continuum background is rejected. 
For the K*~j (K*°j) mode we use the selection criteria 
for the p~-f (p°7) mode. 

We consider the following backgrounds from B decays: 
B — > K*j, other B — » X s j processes, decays with a 7r /?/ 
(B — > /57T°, dJ7r°, prj and wrj), other charmless hadronic 
B decays, and b — > c decay modes. We find the b — > c 
background to be negligible. The B — > K*j background 
can mimic the B — > pj signal if the kaon from the K* is 
misidentified as a pion. To suppress B — ► K*j events we 
calculate Mk-k , where the kaon mass is assigned to one of 
the charged pion candidates, and reject the candidate if 
Mk-k < 0.95 (0.92) GeV/c 2 forthep°7 (p-7) mode. This 
requirement removes 82% (64%) of the K*-f background 
while retaining 63% (87%) of the signal. The decay chain 
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B" -> if* u 7 , -> K s ir°, K s -> tt+tt- has a small 
contribution to B° — > W7 due to the tail of the if* Breit- 
Wigner lineshape. In addition, B — ► if *7 and other £> — > 
X s 7 decays contribute to the background when the p and 
lo candidates are formed from random combinations of 
particles. 

Hadronic decays with a tt° /rj can mimic the signal if a 
photon from the ir° or 77 — > 77 decay is soft and passes 
the 7r°/77 veto. To suppress this background, we reject 
the candidate if | cos Shell > 0.75, 0.70 and 0.80 for the 
p~7, /9°7 and ujj modes, respectively, where the helicity 
angle #hei is the angle between the 7r~ track (normal to 
the uj decay plane) and the B momentum vector in the 
p (u>) rest frame (similarly for the If* 7 modes). Other 
hadronic decays make smaller contributions. 

The reconstruction efficiency for each mode is defined 
as the fraction of the signal remaining after all selec- 
tion criteria are applied, where the signal yield is de- 
termined from a fit to the sum of the signal and con- 
tinuum MC samples using the procedure described be- 
low. The total efficiencies are listed in Table [I] The 
systematic error on the efficiency is the quadratic sum 
of the following contributions, estimated using control 
samples: the uncertainty in the photon detection effi- 
ciency (2.2%) as measured in radiative Bhabha events; 
charged tracking efficiency (1.0% per track) from par- 
tially reconstructed D*+ -> D°n+, D° -> K%-k + -k~ , 
Kg — ► 7r + (7r~); charged pion and kaon identification 
(0.7-1.7% per track) and misidentification (15-17%) from 
D* + — > D a n + , D° — > K~n + ; neutral pion detection 
(4.6%) from 77 decays to 77, 7r + 7r~7r° and 3ir°; TZ-r and 
tt°/j] veto requirements (2.8-5.5%) from B~ — > D°ir~, 
D° -> K-n+ and B° -> D+ir-, D+ -> K-ir+ir+; the 
u) — > 7r + 7T _ 7r° branching fraction (0.8%); and uncertainty 
due to MC statistics (0.5-0.7%). 

We perform an unbinned extended maximum likeli- 
hood fit to candidates satisfying |A.E| < 0.5 GeV and 
Mbc > 5.2 GeV/c 2 , individually and simultaneously for 
the three signal modes. In the latter case we assume 
isospin symmetry, and we also simultaneously fit the two 
B — > if *7 modes. We describe the events in the fit re- 
gion using a sum of functions for the signal, continuum, 
K*j (for the three signal modes only), and other back- 
ground hypotheses. The signal distribution is modeled 
as the product of a Crystal Ball lineshape 0] m 
to reproduce the asymmetric ECL energy response, and 
a Gaussian (another Crystal Ball lineshape) in Mb c for 
the mode without (with) a 7r° in the final state. The 
signal parameters for Mb c and AE are determined from 
separate fits to the B~ —> if *~7 and B° — ► if*°7 sam- 
ples for the modes with and without a neutral pion, re- 
spectively. The branching fraction is the only param- 
eter that is allowed to float for the signal component. 
The continuum background component is modeled as 
the product of a linear function in AE and an AR- 
GUS function in Mbc- The continuum shape pa- 



rameters and normalizations are mode dependent and 
allowed to float. We use the distributions of MC events 
to model the shapes of other background components. 
The size of the if *7 background component in each sig- 
nal mode is constrained using the fit to the if* 7 events 
and the known misidentification probability. Other ra- 
diative and charmless decays are considered as an addi- 
tional background component when we extract the signal 
yield. The levels of the other backgrounds are fixed using 
known branching fractions or upper limits [l4j . We con- 
strain branching fractions in the simultaneous fit using 
the isospin relations HEI B(B -> (p,u))j) = B(B~ -> 
p-j) = 2^f±B(B° -> p°j) = 2^B(B° -> ury) and 

B(B -> If* 7) = B(B- -> K*~i) = T -f^B{B Q -> K*°j), 

where ^ = 1.076 ± 0.008 El- 
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FIG. 2: Projections of the fit results to Mbc (in the re- 
gion -0.10 GeV < AE < 0.08 GeV) and AE (in the region 
5.273 GeV/c 2 < M hc < 5.285 GeV/c 2 ) for the individual and 
simultaneous fits. Curves show the signal (dashed), contin- 
uum (dotted), B — * if* 7 (dot-dashed), other B decay back- 
ground (dot-dot-dashed) components, and the total fit result 
(solid). 



The results of the fits are shown in Fig. and listed in 
Table U The simultaneous fit gives 



b{b^ (p, W ) 7 ) = (i.32 ±S:ii ±8:52) x 



(1) 



where the first and second errors are statistical and sys- 
tematic, respectively. The result is consistent with previ- 
ous results and in agreement with SM predictions 0, 



5 



TABLE I: Yield, significance with (without) systematic un- 
certainty, efficiency, and branching fraction (B) for each mode. 



Mode 


Yield 


Signif. 


Efficiency (%) 


B (10~ 6 ) 


B° -v o>7 
B -» (p,w)7 


8.5 
20.7 
5.7 
36.9 


1.6 (1.6) 

5.2 (5.2) 

2.3 (2.6) 
5.1 (5.4) 


3.86 ± 0.23 
4.30 ±0.28 
2.61 ± 0.21 


n cc +0.42 +0.09 
u -°° -0.36 -0.08 
i 9 c +0.37 +0.07 
-0.33 -0.06 
+0.34 +0.05 
U.OU _ .27 -0.10 
i on +0.34 +0.10 
1 - oz -0.31 -0.09 



5]. The significance of the simultaneous fit is 5.1a, 
where the significance is defined as yj— 21n(£o/X max ), 
and £ ma x (£o) is the value of the likelihood function 
when the signal branching fraction is floated (set to zero). 
Here, the likelihood function from the fit is convolved 
with a Gaussian systematic error function in order to in- 
clude the systematic uncertainty. The invariant 7777(77) 
mass and helicity angle distributions for the events in 
the signal box are consistent with those expected from 
the sum of the signal and background components. The 
fit also gives B(B -> K*j) = (41.1 t\i) x 1(T 6 (statisti- 
cal error only), which is consistent with the world aver- 
age value [l3 |. The individual fit results are in marginal 
agreement with the isospin relation. We test our fitting 
procedure using MC simulation and find no statistically 
significant bias. We perform MC pseudo-experiments 
where events are generated according to the isospin re- 
lation; from the two-dimensional distribution of the de- 
viation between the B~ and averaged B° rates and that 
between the p°7 and ivy rates, we find the probability to 
observe an isospin violation equal to or larger than our 
measurement to be 4.9%. The expected level of isospin 
violation is within ±10% 0. 

The systematic error is estimated by varying each 
of the fixed parameters by ±lcr and then taking the 
quadratic sum of the deviations in the branching frac- 
tion from the nominal value. We note that the ARGUS 
background shape in the fit to the W7 mode is steeper 
than those for the other two modes. Therefore we also 
vary the ARGUS shape parameter for the o>7 mode by 
—2cr and include the deviation in the systematic error. 

The ratio B{B -► (p, u)j)/B(B -> K*j) = 0.032 ± 
0.008(stat.) ± 0.002(syst.), which we obtain from a sep- 
arate fit, can be used to determine |T4d/Vt s |- The fit 
takes into account the correlated systematic errors be- 
tween the signal and K*j modes and thus gives a re- 
duced total error. Using the relation Eal B ^ B ^^ ?hl — 

" """^ B(B— >K*7j 

2 a-"^)/"^) 3 



-C 2 [l ± Ai?], where the form factor 



Vtd 

ratio £ = 0.85 ± 0.10 and the S'C/(3)-breaking correction 
AR = 0.1 ±0.1, we obtain 



\Vu/V ts \ = 0.199^(exp.)™«(theo.). (2) 

We obtain a 95% confidence level interval of 0.142 < 
|Vtd/Vt s | < 0.259 using an ensemble of MC samples in 



which the experimental error is a quadratic sum of the 
asymmetric Gaussian statistical and systematic errors, 
and the theory error is a flat distribution in the given 
range. This result is in agreement with the range favored 
by a fit to the unitarity triangle 01 assuming \Vt s \ — 

\Vcb\- 

In conclusion, we observe the process b —> d'y using 
the B — > p7 and W7 modes. The resulting branching 
fractions are consistent with SM predictions [J, |5j . The 
ratio of the B — > (p, w)j branching fraction to that for 
B — > K*j is used to determine |V t d/Vt s |. 
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